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We perform the two—point diagnostic for the Om{z) fnnction proposed by Sahni etal in 2014 for 
the Starobinsky and Hu & Sawicki models in f{R) gravity. We show that the observed values of 
the Omh^ function can be explained in f{R) models while in LCDM the Omh^ funticon is expected 
to be a redshift independent number. We perform the analysis for some particular values of 
founding a cumulative probability {P{x^ < Xmodei)) P ~ 0.16 or ~ 0.09 for the better cases versus a 
cumulative probability of P ~ 0.98 in the ACDM scenario. We also show that these models present 
a characteristic signature around the interval between z ~ 2 and z ^ A, that could be confronted 
with future observations using the same test. 
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I. INTRODUCTION 

At the present the late time accelerated expansion of 
the Universe is a widely accepted fact supported by dif¬ 
ferent independent observations Recently, several 

modifications of gravity have been proposed in order to 
explain this acceleration, however it is not clear if behind 
these generalizations of the Einstein-Hilbert action there 
is a modification or extension of some physical principle. 
Therefore, there is not a formal framework to guide the 
exploration of these models. 

Every modification of gravity should provide similar 
predictions to those of the ACDM model, nevertheless 
this occurs only at some scales or in certain scenarios, 
e.g., the expansion of the Universe at low 0 or in the 
Solar system tests. These modifications are different from 
each other, so every one of them should provide some 
signature to be distinguished from the others. Because 
of the lack of a formal guidance, the principal line of 
exploration for any model of modified gravity is to impose 
observational constraints on their parameters relying on 
the high accuracy of the current data. 

The /(i?) gravity constitute a natural extension of 
General Relativity in the sense that almost all of its de¬ 
sirable properties remain intact. One of the changes in¬ 
troduced by the dependence of the gravitational action 
on a general function of the Ricci scalar, f{R), is that 
we increase the number of degrees of freedom, resulting 
in a different set of field equations and consequently in 
modified dynamics for the gravitational systems. [3-@| 

Among the components of these modified gravities 
which can be tested are the equation of state (EOS) of 
the dark energy component and the growth factor. The 
values used in performing these tests can be obtained, in 
an indirect way, from astronomical measurements. These 
analyses depend on the value of the matter density today 
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for both dark and barionic components and even a 
small change can produce significant differences in the 
behavior of these models. It is worth mentioning that 
for f{R) gravity an additional problem emerges, that is 
the necessity of a unique definition for the EOS of the 
geometric dark energy [l^ . 

In these context, Sahni et al. EH proposed a new test 
in order to distinguish between acceleration produced by 
a cosmological constant term and that coming from a 
modified gravity model. Now that results at different 
redshift for the Hubble parameter H{z)/H{z = 0) are 
available from the observations of Barionic Acoustic Os¬ 
cillations (BAO), it is possible to implement this test for 
AGDM as well as for modified theories of gravity. 

Among the more successful models in f{R) we find 
that of Starobinsky [l^ and the proposed by Hu & Saw¬ 
icki [ 1 ^. These models have proved to be viable at low 
redshift z ~ 1. In the present work we analyze them 
to verify their applicability at a higher redshift z Ri 2.34 
which is is within the scale of validity for tests using the 
BAO data. 

In the next section we will present the approach that 
we are using [I^ in order to integrate the modified Fried¬ 
man equations and the /(i?) models we will explore. In 
section HI the Omh{z) in introduced as well as the two- 
ponit relation Omh^^Zi] Zj) that is used as a model inde¬ 
pendent test in this paper. Section IV shows the results 
for the Ohm diagnostic for the f{R) models and they are 
compared with the AGDM model. Finally, in section V 
we present teh conclusion of this work. 


II. f{R) COSMOLOGY 

As we have mentioned, f{R) theories of gravity are the 
most straightforward way to extend the Hilbert-Einstein 
action. The dependence of the Ricci scalar is a general 
function which will be defined in order to reproduce ob- 



servations, the action is given by 


• Starobinsky (l^ : 


S[gab, 4’] 



•S^inatt 


[5ah,'0] , 


( 1 ) 


where G = 1, c = 1 and k = Stt. The term f{R) 
is an arbitrary smooth function of the Ricci scalar and 
Smatt [dab, 4’] is ttie usual action for matter. 

Varying the above action with respect to g°‘^ we obtain 
the modified field equations 

fnRab — -^fgab — (VaVfe — gab^) fu = bcTab , (2) 


where Jr = OrJ, □ = g’^'^S/a^b and Tab is the energy- 
momentum tensor for matter. The set of equations can 
be re-written in the following way 


fuGab - fRRVaVbR " fRR r{V aR){VbR) 


+ gab 


Wr -f) + IrrOR + fRRRi^Rf 


— RTab ; 

(3) 


where Gab = Rab — gabR/‘2. is the Einstein tensor and 
(Vi?)2 := g-\VaR){VbR). 

The trace of eq. dU yields a second order equation for 
the Ricci scalar 


□R = 


1 

3/fifl 


- 3/flflfl(VR)2 + 2/ - Rfn 


where T := T“. Finally, using (|1]) in ([3]) we find 


(4) 


Gab — - 7 — 
JR L 
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fRR^a^bR + /flflii(VaR)(Vf,i?) 


{rJr 


+ / + 2kT ) + uTab 


(5) 


In this work we will consider an homogeneous, isotropic 
universe described by the Friedman-Robertson-Walker 
metric 


ds^ = —dt^ + a^{t) 


dr^ 


1 — kr'^ 


T (d6»2 -p sin^ Odtp^) 


( 6 ) 

as usual, we will assume k = 0. The energy momentum 
tensor (EMT) is that for a fluid composed by baryons, 
dark matter and radiation. Under these assumptions 
Eqs. dH) and dSI read 


R= -3HR- 


1 


^Irr - 


SfRRRR^ + 2f- /rR + kT (7) 


1 


H^ = -^ IrrHR--{RJr- f) 
JRR L 6 


H= -H‘^- 


H = d/c 


k-T* 

(9) 

( 10 ) 


where ' = d/dt. We will consider the following f{R) 
models: 


f{R) 


R -\- XRs 



with g = 2, A = 1 and Rs = 4.17 Hq. 
• Hu-Sawicki (l5j |: 


(11) 


/(R) = R 



( 12 ) 


where the parameters are n = 4, ci ~ 1.25 x 10 ^ 
C 2 « 6.56 X 10“^ and Rhs ~ 0.24Rg. 


These two models are currently the most successful 
among the proposed /(R) modified gravity descriptions 
and both pass the Solar System tests. 

The numerical integration will be performed by using 
equations © and dH) together with the standard conse- 
vation equation pi + 3H (pi + Pi) = 0 with i = 1,..., 3 
for each fluid component, namely baryons, dark matter 
and radiation. The Hamiltonian constriction eq. m is 
used in order to check the error in the numerical code 
which is acceptable (^ 10“^^). 


III. THE Omh DIAGNOSTIC. 


Sahni et al. proposed [ll| a test to distinguish ACDM 
from modified gravity models or some other mechanism 
to describe the late acceleration of the Universe. The test 
is based on the function Om defined as 


Om(x) = 


h'^{x) — I 
x^ — 1 


(13) 


with X = 1 + z and h(x) = H(x)/Hq. It is worth men¬ 
tioning that the Om function is constant for an accel¬ 
erated expansion described by means of a Cosmologi¬ 
cal Constant term instead of a modified gravity model. 
In fact, in such case we have that Om{x) = thus 
Om{x) — = 0 and any deviation from zero would dis¬ 

card ACDM as a model for the expansion of the Universe. 
In a modihed gravity description Om{x) is not constant 
but evolves with 2 , showing a particular behavior for each 
theory. 

A remarkable characteristic of this test is that it de¬ 
pends only on the values of H{z) which are determined 
by observations. 

Shaheloo et al. 


I7| proposed a diagnostic by using 


the Om function at two different points. This way we 
can take observations about the determination of H{z) 
at several redshifts and then we can compute the two 
point relation given by 
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(14) 


0.32 


Om{z2;zi) 


fe^(z2) - h’^jzi) 

(1 + Z2)3-(1 + Zi)3' 


Later, Sahni et al. in [T^ proposed a small modi¬ 
fication to this test which allows to compare different 
theoretical predictions using the two-point diagnostic of 
Shafieloo et al. They multiply Eq. (imi by with 
h = HQ/lQQkm/sec/Mpc obtaining 


Omh'^{zi; zj) 


h'^izi) - li^izj) 

(l + Z.)3-(l + ^,)3- 


(15) 


Written in this forms, this little change have the ad¬ 
vantage that we can use that under the ACDM model 
riMh"^ = 0.1426 ± 0.0025 from Planck XVI 2013 [H. 
Sahni et al. [l^ show that, by using this test, ob¬ 
servations suggest that the value of Omh^ is not con¬ 
stant. Taking zi = 0, ^2 = 0.57 and z^ = 2.34 with 
H{zi) = 7Q.%±i.2kmlseclMpc [Hi^, H{z 2 ) = 92.4± 
A.bkm/sec/Mpc [2l| and H{z^) = 222 ± IkmjsecjMpc 
[T^ the values for the Omh"^ two-point relation (eq. [TS] 1 
reported by Sahni are 


Omh'^{zi-,Z2) = 0.124 ± 0.045, 

Om/i2(zi;z3) = 0.122 ±0.01, (16) 

Omh‘^{z2-, z^) = 0.122 ± 0.012, 

while for ACDM the value is Omh^ = 0.1426 with 
a cumulative probability P{ACDM) = 0.98, i.e. the 
probability to find a sample with a < Xmodei i^ix^ — 
xLdez))> so the p-value is given by 1 - ^(x^ < xLdez)- 
We have considered just this value for the ACDM model 
as this is the value obtained in the best fit given by the 
Planck team (T9j |. 

In the next section we apply this diagnostic to two 
f{R) models. We will present the Om{z) function, the 
two-point relation Omh^{zi; Zj) will be also computed in 
bothj cases. f{R) results will be compared with ACDM 
in terms of the cumulative probability. 


IV. f{R) DIAGNOSTIC 

We have intruduced by now the way we will integrate 
the field equations in cosmology and also the test we will 
use as a diagnostic for two f{R) models. To perform 
the analysis we proceeded as follows. We integrate the 
differential equations o and under the Ricci scalar 
approach fixing the initial conditions at some point in the 
past where it is safe to assume domination of matter and 
normalizing in order to fix the values of (for details 
about cosmological integration in the context of the Ricci 
scalar approach see ). It is important to mention that 
we are taking into account a radiation component in the 
EMT. Nevertheless, this has no a significant effect in the 



FIG. 1. (colour online). Om{z) function for the Starobin- 
sky model assuming (top to bottom) = 0.286 (red), 0.26 
(blue), 0.25 (orange), 0.24 (gray) and 0.23 (black). 


results when integrating in the periods we are consider¬ 
ing, i.e. z < 5. 

The value of assumed in the theoretical predic¬ 
tion, is crucial in order to have a better agreement with 
the values of the two-point relation given in (1161) . In any 
case it is clear that, even if such changes are considered, 
the ACDM model can only give constant values for the 
Om{z) function, and consequently it gives constant val¬ 
ues for the two-point relation (IT4l) . 

Here we explore different values of for each f{R) 
model: = 0.286 (which corresponds to the value ex¬ 

pected under the ACDM model) and also = 0.23, 
0.24, 0.25 and 0.26. In Eig. 1 we show the evolution of 
Om{z) given by Eq. (ITSl) as a function of the redshift 
for the Starobinsky model assuming different values for 
In Eig. 2 we present the same function for the Hu 
& Sawicki model. The two models have a peculiar be¬ 
havior that is worth mentioning. For a redshift around 
the value z ^ 2 the function Om(z) reaches a minimum 
and from there it takes an stable or nearly stable value 
(see Fig. |31). This behavior shows a particular predic¬ 
tion from these models of f{R) gravity that could be 
confronted with observations via the two-point relation. 

Therefore, for each value we compute the two 
point relation Omh?{z\]Z 2 ) given by Eq. (ITSl) . The val¬ 
ues obtained for the Starobinsky model are presented 
in Table I. In Table H we show the values of the two 
point-relation Omli^{zi; Z 2 ) obtained for the Hu & Saw¬ 
icki model. As we mentioned, in the ACDM model, this 
value is constant Omh? = 0.1426. 

The cumulative probability for each model {P{Model)) 
has been computed taking the x^ values in tables I and 
H. The values are incorporated in column 3 of Tables I 
and H. We can observe the best case for the Starobinsky 
model is for = 0.24 with a cumulative probability 
Pif{R)st) = 0.16. For the Hu & Sawicki model the 
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{Zi,Zj) 

Omh‘^{zi, Zj) 

Xst 

P{fiR)st) 


[Z1,Z2) 

0.145 



0.286 

(zi, 23 ) 

0.139 

4.978 

0.97 


( 22 , 2 : 3 ) 

0.138 




( 21 , 22 ) 

0.137 



0.26 

( 21 , 23 ) 

0.130 

1.179 

0.72 


( 22 , 23 ) 

0.130 




( 21 , 22 ) 

0.134 



0.25 

( 21 , 23 ) 

0.127 

0.399 

0.47 


( 22 , 23 ) 

0.126 




( 21 , 22 ) 

0.131 



0.24 

( 21 , 23 ) 

0.123 

0.041 

0.16 


( 22 , 23 ) 

0.123 




( 21 , 22 ) 

0.128 



0.23 

( 21 , 23 ) 

0.120 

0.136 

0.28 


( 22 , 23 ) 

0.119 




TABLE I. Values for the two points two-point relation Omh^{zi, Z 2 ) for the Starobinsky f{R) model assuming different values 
of Column 1 shows the value of the prior for the value. In column 2 we have the two points in 2 : where Omh^{zi, Zj) 
is calculated (21 = 0, 22 = 0.57 and 23 = 2.34). In column 3 is the value of Omh. In column 4 is the value of x%f In column 5 
is the cumulative probability. 



FIG. 2. (colour online). Om{z) function for the Hu & Saw- 
icki model assuming (top to bottom) Q%[ = 0.286 (red), 0.26 
(blue), 0.25 (orange), 0.24 (gray) and 0.23 (black). 
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FIG. 3. Evolution of the Om{z) function. Starobinsky and 
Hu & Sawicki models, top and bottom respectively, assuming 
0,% = 0.25 in both cases. 


best case is when we have = 0.25 with a cumulative 
probability P{f{R)H-s) = 0.09. 


V. CONCLUSION 


In the ACDM case = 7.361 and the cumulative 

probability P{ACDM) = 0.98, these values are not listed In this work we performed the Omh^ test for two of 
in Tables I and II. the most successful f{R) models taking several values 
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‘‘m 

{Zi,Zj) 

Omh‘^{zi, Zj) 

Xh-s 

P{f{R)H-s) 


( 21 , 22 ) 

0.138 



0.286 

( 21 , 23 ) 

0.136 

3.144 

0.92 


( 22 , 23 ) 

0.135 




( 21 , 22 ) 

0.129 



0.26 

( 21 , 23 ) 

0.126 

0.328 

0.43 


( 22 , 23 ) 

0.126 




( 21 , 22 ) 

0.126 



0.25 

( 21 , 23 ) 

0.123 

0.013 

0.09 


( 22 , 23 ) 

0.123 




( 21 , 22 ) 

0.122 



0.24 

( 21 , 23 ) 

0.119 

0.137 

0.29 


( 22 , 23 ) 

0.119 




( 21 , 22 ) 

0.118 



0.23 

( 21 , 23 ) 

0.116 

0.750 

0.61 


( 22 , 23 ) 

0.115 




TABLE II. Values for the two points two-point relation Omh^{zi, 22 ) for the Hu-Sawicki f{R) model assuming different values 
of Column 1 shows the value of the prior for the value. In column 2 we have the two points in z where Omh^{zi, zj) 
is calculated (21 = 0, 22 = 0.57 and 23 = 2.34). In column 3 is the value of Omh. In column 4 is the value of x|t- In column 5 
is the cumulative probability. 


of the current density 17^. The results obtained show 
that the two considered f{R) models have, in general, 
the behaivour that is expected from observations. In the 
case of ACDM this behaivour can not be present because 
Omh^ is expected to be a redshift independent number. 

In the /(i?) models we use in this paper, the evolution 
of the Om{z) function is appropiate to have the values 
expected under the two-point relation, as we can observe 
for all the elections of 17^, it is found that the cumu¬ 
lative probability is better than in the ACDM case. In 
particular, for the Starobinsky model with = 0.24 
the cumulative probability P{f{R)st) = 0.16 and for 
the Hu-Sawicki case we find P{f{R)H-s) = 0.09 tak¬ 
ing = 0.25. This both results are suitable values for 
this type of statistical test. 

The asymptotic behavior of Omh^ in these f{R) mod¬ 
els is necessary to reach the values from Omh‘^{zi, Z3) 


to Omh?{z 2 , Z 3 ) which remain (almost) constant. This 
behavior is not possible using the ACDM model. It is re¬ 
markable that the evolution of the Omh? function from 
2 ~ 2 to 2 ~ 4 could be tested by using future observa¬ 
tions. Such evolution does not decrease in a monotonic 
way, it presents a change of sign of the Omh‘^{zi, Zj) and 
this signature could be used in order to rule out or give 
support to these models. 
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